Polymer translocation out of confined environments 
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We consider the dynamics of polymer translocation out of confined environments. Analytic scaling 
arguments lead to the prediction that the translocation time scales like r ~ j\[I 3 + v 2d jj1+(i-"2d)/ !/ f or 
translocation out of a planar confinement between two walls with separation R into a 3D environ- 
ment, and r ~ N^ +1 R for translocation out of two strips with separation R into a 2D environment. 
Here, N is the chain length, v and i*2d are the Flory exponents in 3D and 2D, and (5 is the scal- 
ing exponent of translocation velocity with TV, whose value for the present choice of parameters 
is j3 ~ 0.8 based on Langevin dynamics simulations. These scaling exponents improve on earlier 
predictions. 

PACS numbers: 87.15.A-, 87.15.H- 



I. INTRODUCTION 

The transport of a polymer through a nanopore has 
received increasing attention due to its importance in bi- 
ological systems |l| , such as gene swapping through bac- 
terial pili, m-RNA transport through nuclear pore com- 
plexes, and injection of DNA from a virus head into the 
host cell, etc. Moreover, translocation processes have 
various potentially revolutionary technological applica- 
tions j such as rapid DNA sequencing, gene therapy, 
and controlled drug delivery. 

In addition to its biological relevance, the transport 
dynamics of polymers through nanoporcs is of funda- 
mental interest in physics and chemistry. There ex- 
ists a flurr y o f experimental [ 41, [ a , [fj, hfl and theoreti- 
cal i, Jl M EL 51 H HliH EE EE HE 

El HE IHTIITIITMTIETII, studies devoted to 
this subject. The passage of a flexible chain through 
a nanopore involves a large entropic barrier, thus poly- 
mer translocation needs driving forces, which can be pro- 
vided by an external applied electric field in the pore 
@, II IE IE EE HE HI [2|j2|La milling force exerted 

29[, binding particles 



on the end of a 
(chaperones) [Ml |1 
polymer fj\, IE Illl.ll 

Among these mechanisms, less attention has been paid 
to confinement-driven translocation. In particular its dy- 
namics remains unclear, despite its importance to bio- 
logical processes including viral ejection, drug delivery, 
controlled release from a nanochannel, etc. We here in- 
vestigate the generic behavior of polymer release through 
a small pore from a confined environment (Fig. [1]). 

Based on Kantor and Kardar's scaling arguments for 
unhindered motion of the chain [ItJ , Cacciuto and Lui- 
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jten [lj| suggested that r ~ N 1+V R}/ V for planar con- 
finement shown in Fig. [1] where v is the Flory expo- 
nent [3l], HE HH (y = 0.588 in 3D) and R the separation 
between the walls. Panja et al. [2(| considered transloca- 
tion out of planar confinement as two-dimensional version 
of translocation through a nanopore with an electric field 
applied in the pore, and suggested r ~ jV 2ly2D (1/2 d = 3/4 
in 2D), in contrast to above scaling exponent 1 + v. To 
solve this contradiction, in this work we investigate the 
dynamics of polymer translocation out of planar confine- 
ment (3D) and out of two strips (2D) using both analytic 
scaling arguments and Langevin dynamics simulations. 
As in the above mentioned theories, we also ignore the 
hydrodynamic interactions here. Our scaling arguments 
include geometric effects that have been left out in previ- 
ous studies, leading to a new set of scaling exponents of 
translocation time r with regard to polymer length N as 
well as the size of the confinement space R. The theoret- 
ical findings for R dependence of r are fully supported by 
numerical simulation results both in 2D and 3D. We con- 
sider a geometry as shown in Fig. [TJ where two walls with 
separation R are formed by stationary particles within a 
distance a from each other. 



II. SCALING ARGUMENT 

A number of recent theories [E EE EE E3 have been 
developed for the dynamics of polymer translocation. Of 
them, Kantor and Kardar [l7j provided lower bounds for 
both the translocation time for pulling the polymer by 
the end and driving the polymer by a chemical poten- 
tial difference applied across the membrane. Essentially, 
the lower bound is the time for the unimpeded polymer 
moving a distance of order of its size. For driven translo- 
cation with a chemical potential difference applied across 
the membrane, there is a force of A/i/a applied to the 
few monomers inside the pore. As a consequence, Kan- 
tor and Kardar [13 argued that its center of mass should 
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FIG. 1: A schematic representation of the system in 2D. The 
simulations are carried out in both 2D and 3D. For a planar 
confinement (3D), two plates are separated by a distance R. 
One wall has a pore of length L — a and diameter W = 3cr. 
For polymers confined between two strips (2D) the pore width 
is 3(J. 



move with a velocity v ~ A/jl/N. Thus the lower bound 
for the translocation time of an unhindered polymer is 
the time to move through a distance R g (radius of gyra- 
tion of the polymer in an unconfined state), and scales 
as r - R g /v ~ iV 1+1 7(A/i). 

Now, let us focus on the translocation out of con- 
fined environments. According to the blob picture [3l| . 
a chain confined between two parallel plates with sepa- 
ration a <C R <C R g will form a 2D self-avoiding walk 
consisting of rib blobs of size R. Each blob contains 
g = (R/a) 1 '" monomers and the number of blobs is 
n b = N/g = N(R/ay 1 ^. The free energy cost of the 
confined chain in units of ksT, AF, is simply the number 
of the blobs. Thus the chemical potential difference that 
provides a driving force for the translocation process is 
given by 



A// = AF/N ~ R 



-l/u 



(1) 



The remaining ingredients to complete the scaling argu- 
ment are the length scale Lq through which the polymer 
has to move during the entire translocation process and 
the average translocation velocity v. As a result, the 
translocation time can be estimated as t ~ Lq/v. Cac- 
ciuto et al. [l9j chose the radius of gyration R g in the 
unconfined state for Lq and v ~ A/x/N and obtained the 
scaling result r - N^R 1 ^ ~ Ar 1 - 588 ^ 1 - 70 , while Panja 
et al. [2(| criticized the scaling exponent 1 + v for N 
dependence of r and considered the translocation out of 
the planar confinement as the 2D version of translocation 
driven by a chemical potential difference applied across 
the membrane. As a result, they obtained the exponent 
2^2D = 1-50 for the N dependence in a 3D system. 

Here, we argue that the correct scaling results can only 
be obtained by properly identifying the length scale Lq 
and the scaling of the translocation velocity v. Due to the 
highly non-equilibrium nature of the driven translocation 
problem, the scaling of the average translocation velocity 
v with respect to the chain length N is of the from v ~ 
N@ , where the exponent f3 < 1 [30| . For the present 
choice of parameters, we find that j3 s=s 0.8 as will be 



demonstrated below. For planar confinement (3D), the 
blob pict ure p redicts the longitudinal size of the polymer 
to be 0,11 1! 
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(2) 



with V2D being the Flory exponent in 2D. For successful 
translocation out of planar confinement, the polymer has 
to move a distance Lq ~ Ru . Therefore, the translocation 
time can be estimated as 

Rn 



V 



(3) 



Numerically, Eq. (|3|) yields a scaling result of r ~ R 1A3 
for confinement driven translocation in 3D, which is dif- 
ferent from the prediction r ~ R 170 obtained by Cacci- 
uto et al. [l9l |. In fact, their numerical results based on 
Monte Carlo simulations show r ~ ^i-54±o.io^ w j 1 j c j 1 j s 
close to our scaling prediction in Eq. ([3"|). 

Next, we also consider translocation out of confinement 
into a 2D environment which has not been addressed pre- 
viously 0, HH, H, H3] • For a polymer confined between 
two strips embedded in 2D, the blob picture predicts the 
longitudinal size of the chain to be [3l|, HH 



R 



(4) 



In this case the free energy excess in Eq. ([T]) is valid if v 
is replaced by V2D- Thus, the translocation time scales 
as 

7?ii N Rl— V^a-D 

T ~ til „ ^ „ N P+i R; (5) 



v 



showing a linear dependence on R. Both predictions in 
Eqs. j3]) and ([5|) for R dependence are confirmed by sim- 
ulations below. 



III. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard-Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction between 
monomers is modeled by a short range repulsive LJ po- 
tential: U LJ {r) = 4e[(f ) 12 - (^) 6 ] + e for r < 2 1 / 6 ct and 
for r > 2 1 / 6 cr. Here, a is the diameter of a monomer, 
and e is the potential depth. The connectivity between 
neighboring monomers is modeled as a FENE spring with 

U fene(t) = — ^fci?o m (l ~ 7,2 /^o)' where r is the dis- 
tance between consecutive monomers, k the spring con- 
stant, and Ro the maximum allowed separation between 
connected monomers. Between all monomer-wall parti- 
cle pairs, there exists the same short range repulsive LJ 
interaction as described above. 

In the Langevin dynamics simulation, each monomer is 
subjected to conservative, frictional, and random forces, 
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FIG. 2: Translocation time r as a function of the chain length 
N for different R in 3D. The insert shows r as a function of 
R for chain length N = 128 in 3D. 

respectively, with [34[ mfj = —V{Ulj + Ufene) — + 
F^, where m is the monomer's mass, £ is the friction co- 
efficient and Ff is the random force which satisfies the 
fluctuation-dissipation theorem. In the present work, we 
use the LJ parameters e and a and the monomer mass 
in to fix the energy, length and mass scales respectively. 
Time scale is then given by t^j = [ma 2 j e) 1 ! 2 . The di- 
mcnsionless parameters in our simulations are Rq = 2, 
k = 7, £ = 0.7 and fc^T = 1.2. The Langevin equation is 
integrated in time by a method described by Ermak and 
Buckholz [HI in both 3D and 2D. To create the initial 
configuration, the first monomer of the chain is placed in 
the entrance of the pore, while the remaining monomers 
arc initially squeezed into the space between two plates 
(3D) and the space between two strips (2D) under ther- 
mal collisions described by the Langevin thermostat to 
obtain an equilibrium configuration. Typically, we aver- 
age our data over 1000 independent runs. 

According to the definition of the translocation time, 
at the completion of the translocation process, the chain 
has moved a distance of Ru along the direction per- 
pendicular to the axis of the pore. We now define the 
translocation velocity as v — (Ru)/t. Using the defini- 
tion v = (i?y./r,) for the translocation velocity, we ob- 
served similar results. Here iJy and Tj denote the values 
of -Ru and t for each successful run. 



IV. NUMERICAL RESULTS 

For translocation out of a planar confinement, the re- 
sults are shown in Fig. [21 The insert of Fig. [2] shows 
the R dependence of the translocation time. The scal- 
ing exponent is 1.42 ± 0.02, which is in good agreement 
with the exponent 1.43 from our scaling prediction in 



Eq. ©. For r with N for R = 3.5, 2.0 and 1.5, we 
get the scaling exponents of 1.46 ± 0.02, 1.48 ± 0.01, and 
1.49 ± 0.03, respectively With decreasing R, the expo- 
nent slightly increases. As to translocation velocity, such 
as for R = 1.5 we get (3 = 0.77 ± 0.01. According to Eq. 
©, (3 + v w = 1.52, which is very close to 1.49. These 
results demonstrate that the scaling arguments for unim- 
peded translocation provide an accurate estimate for the 
translocation time. Although the scaling exponent for N 
dependence is quite close to the value 1.55 obtained by 
Cacciuto et al. [191 ]. the physical origin is quite different. 

The dynamics of a single segment passing through the 
pore during translocation is an important quantity con- 
siderably affected by different driving mechanisms. The 
noncquilibrium nature of translocation has a significant 
effect on it. We have numerically calculated the waiting 
times for all monomers in a chain of length N. We de- 
fine the waiting time of monomer s as the average time 
between the events that monomer s and monomer s + 1 
exit the pore. In our previous work [25L [26j studying 
translocation driven by a voltage across the pore with- 
out confinement on either side, we found that for short 
polymers, such as N — 100, the monomers in the middle 
of the polymer need the longest time to translocate and 
the distribution is close to symmetric with respect to the 
middle monomer [26| . Fig. [3] shows the waiting time dis- 
tributions for translocation out of planar confinements 
with R = 2.0 and 3.5 for N = 64. Compared with the 
unconfined potential driven case, the waiting time distri- 
bution is quite different. The waiting times increase until 
s ~ 20, and then almost saturate. This can be qualita- 
tively understood as follows. The initial rise of the wait- 
ing time has the same origin as the potential driven case, 
with A/i given by Eq. ([T]) playing the role of the applied 
voltage. Here, the driving force has two components, 
one from the cost of the free energy due to blobs and the 
second due to entropic wall repulsion for the translocated 
beads. These are balanced by frictional force in the pore. 
The reduction in the overall driving causes the initial rise 
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FIG. 4: Translocation time r as a function of the chain length 
N for R — 4.5 in 2D. The insert shows r as a function of R 
for chain length N = 128 in 2D. 



linearly with increasing R with the scaling exponent 
1.04 ± 0.01, see the insert of Fig. QJ which is in excel- 
lent agreement with the prediction in Eq. ([5]). More- 
over, for r function of N, we get scaling exponent 
of 1.82 ± 0.01 for R = 4.5. For translocation velocity, we 
get /3 = 0.79 ± 0.02. According to Eq. © , /? + 1 = 1.79, 
which is very close to the scaling exponent 1.82. 

Fig. [5]shows the waiting time distributions for translo- 
cation out of two strips with R = 4.5 for N = 16, 32, 64 
and 96. For N > 64, the waiting time increases with s 
until the maximum at s ~ 20, and then almost linearly 
decreases with s. For all chain lengths, the monomers at 
the end of the chain from s = N — 4 to N — 1 still need a 
little longer time to translocate due to very weak confine- 
ment. The observed behavior is due to the interplay of 
decreasing Afj, and increasing pulling entropic force. The 
balance of these two factors is differs from 2D to 3D. 



V. CONCLUSIONS 



R=4.5 
□ N=16 
o N=32 
a-N=64 
v N=96 




In this work, we have considered the dynamics of poly- 
mer translocation out of confined environments based 
on both scaling arguments and Langevin dynamics sim- 
ulations. Analytic scaling arguments lead to the pre- 
diction that the translocation time scales like t ~ 
jy/3+f2D#i+(i-i/2D)/V f or translocation out of a planar 
confinement between two walls with separation R into 
a 3D environment, and r ~ N@ +1 R for translocation out 
of two strips with separation R into a 2D environment. 
Here, N is the chain length, v and vid are the Flory ex- 
ponents in 3D and 2D, and P is the scaling exponent of 
translocation velocity with TV, whose value for the present 
choice of parameters is P « 0.8 based on Langevin dy- 
namics simulations. These scaling exponents improve on 
the previously reported results [19|, [20j ■ 



FIG. 5: Waiting time distribution for 2D simulations. 
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